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LIMIT DISTRIBUTIONS OF GAUSSIAN BLOCK ENSEMBLES 


ROMUALD LENCZEWSKI 


Abstract. It has been shown by Voiculescu that important classes of square inde¬ 
pendent random matrices are asymptotically free, where freeness is a noncommutative 
analog of classical independence. Recently, we introduced the concept of matricial 
freeness, which is similar to freeness in free probability, but it also has some matri¬ 
cial features. Using this new concept of noncommutative independence, we described 
the asymptotics of blocks and symmetric blocks of certain classes of independent ran¬ 
dom matrices. In this paper, we present the main results obtained in this framework, 
concentrating on the ensembles of blocks of Gaussian random matrices. 


1. Introduction 

Our main objective is to describe asymptotic joint distributions of rectangular blocks 
of independent random matrices, called random blocks, under the expectation of nor¬ 
malized partial traces. For that purpose, we use a new concept of noncommutative 
independence called matricial freeness and associated arrays of operators which give 
Hilbert space realizations of these distributions. The well-known connection between 
free probability and the asymptotics of independent square random matrices under the 
expectation of normalized trace can also be reproduced in this framework. 

The most fundamental results of this nature were obtained by Voiculescu [16], who 
showed that certain ensembles of independent n x n random matrices {Y{u,n) : u e 
Zl} were asymptotically free under the expectation of the normalized complete trace, 
r(n) = E (X) Tr(n), where Tr(n) stands for the trace divided by n. In particular, if 
the entries of Y{u,n) are i.i.d. complex (0, l/n)-Gaussian random variables, we can 
symbolically write 

lim Y{u, n) = r]{u) 

n^GO 

where {r]{u) : u e Zl} is the standard free circular system of operators and convergence 
is understood in the sense of mixed moments under T{n). The operators rj{u) live in 
the free Fock space and have the standard circular distribution (uniform distribution 
on the unit disc in the complex plane) in the vacuum state. A similar result holds for 
Hermitian Gaussian random matrices whose limit joint distributions are described by 
mixed moments of free Gaussian operators with semicircle distributions. 

Gaussian random matrices studied by Voiculescu had i.i.d. entries, except that in 
the case of Hermitian ensembles it holds that Yij{u,n) = Yj^i{u,n). If the entries are 
independent but not identically distributed, standard free probability may not suffice 
to describe asymptotic distributions of these matrices. One of the approaches is to 
employ the much more general scheme of freeness with amalgamation as in the papers 
of Shlyakhtenko [14] and Benaych-Georges [3]. Recently, we have developed another 
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approach [7,8] in the case when the entries are independent and block-identically dis¬ 
tributed {i.b.i.d.). This approach is based on the concept of matricial freeness [6], 
which can be viewed as a matricial generalization of freeness, lying somewhere between 
freeness and freeness with amalgamation. 

In this context, we decompose matrices Y{u,n) with i.b.i.d. entries into rectangular 
blocks {Sp^q{u, n) : 1 ^ p,q ^ r}, or symmetric blocks {Tp^q{u, n) : 1 < p < g ^ r}. The 
block analog of the above Voiculescu’s result for the ensemble of independent Gaussian 
matrices can then be written in the form 

lim Sp^q{u, n) = Cp,q{u) 
n^co 

where Qp^q{u) are certain operators living in the matricially free Foek spaee, a matricial 
analog of the free Fock space, into which ri{u) decompose. The convergence is in the 
sense of mixed *-moments under the expectation of normalized partial traces over the 
subsets of basis vectors dehned by the block decomposition, denoted Tq{n) = EoTrg(?7.), 
where 1 < g < r. In turn, the corresponding moments of operators are computed with 
respect to the family of vacuum states. A similar theorem holds for symmetric blocks 
of Hermitian and non-Hermitian Gaussian random matrices. 

This block model encodes more degrees of freedom than the usual framework of 
square random matrices with i.i.d. Gaussian entries. Apart from block variances, the 
most important parameters are the asymptotie dimensions of blocks, namely 

dq = lim —, 

n^oo n 

where n = ni -t- ... -1- n,. is the n-dependent partition of n dehned by the block decom¬ 
position of Y{u,n). We assume that these limits exist, but it is possible that some of 
them vanish. This leads to three types of blocks or symmetric blocks: balaneed (if both 
their asymptotic dimensions are positive), unbalaneed (if one asymptotic dimension is 
positive and the other one vanishes) and evaneseent (if both asymptotic dimensions 
vanish). For details, see [8]. 

Using all types of blocks, we can construct random matrix models for other notions 
of independence, such as monotone independence, boolean independence, c-freeness 
(conditional freeness) and s-freeness (freeness with subordination). Moreover, the limit 
moments can be written in a quite explicit form as polynomials in di,...,dr, with 
block variances as additional parameters. In some cases, this enables us to construct 
simple random matrix models for certain families of probability measures, like the 
model for free Meixner laws [1,4] constructed in [9]. At the same time, studying such 
polynomials can lead to new results in free probability. For instance, a study of the 
limit distributions of products of independent rectangular Gaussian random matrices 
produced polynomials which can be viewed as multivariate Fuss-Narayana polynomials 
and, moreover, turned out to be the moments of the free multiplicative convolution of 
Marchenko-Pastur distributions with arbitrary shape parameters [10]. 

In order to make this paper easy to follow, we avoid technical details and only sketch 
some proofs. For complete proofs, we refer the reader to [7,8,9,10]. 

2. Matricial operator systems 

In the framework of free probability, an important role is played by semicircular 
and circular systems of operators. Such systems were introduced by Voiculescu, who 
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used them to construct the generators of free group factors in order to prove certain 
isomorphisms between them [16]. 

These generators play also an important role in the framework of matricial freeness. 
For the dehnition of that concept, which can be viewed as a matricial generalization of 
freeness involving a family of states instead of one state, we refer the reader to [6]. In this 
paper, we prefer to restrict our attention to important examples of arrays of operators 
[7,8] which are matricially free with respect to a family of states associated with vacuum 
vectors in the ‘matricial’ analog of the free Fock space, called the matricially free Fock 
space. 

Definition 2.1. Let {'Hp^q{u)), u e Zl, he a family of r x r arrays of Hilbert spaces. 
The matricially free Fock space is the Hilbert space direct sum of the form 

q=l 

where 

00 

Afq = CHq 0 0 0 'Hpi,p2{Ul) <S) ■ ■ ■ 

m=l 

ui,...,Un€‘U 

with [r] := {l,2,...,r} and {Hi,..., being the set of unit (vacuum) vectors, 
equipped with the canonical inner product. Let {Ti,...,^^} be the corresponding 
states replacing the single vacuum state in the free Fock space. 

In fact, for most purposes, it suffices to take each Hilbert space to be one-dimensional, 
namely 'Hp^q{u) = Ccp^g^u), where the vectors Cp^qiu) form an ortho normal basis. Arrays 
of certain partial isometries living in this Fock space which remind free creation oper¬ 
ators i.{u) living in the free Fock space serve as generators of certain Toeplitz-Cuntz- 
Krieger algebras [5]. These partial isometries, when multiplied by positive scalars, 
become matricial analogs of free creation operators. 

Definition 2.2. Let B{u) = {bp^q{u)) be an array of positive real numbers for any ue U. 
We associate with each such array the array of matricially free creation operators whose 
non-trivial action onto the basis vectors is 

Pp,q{u)klq = \J hp^q{u)ep^q{u) 

Ipp,q{u){eq^t{,s)) = hp^q{u){ep^q{u) ® eq^t{s)) 

Pp,q{u){eq^t{s) ® w) = hp^q{u){ep^q{u) ® eq^t{,s) ® w) 

for any p,q,t e [r] and u,s e U, where eq^t{s) 0 tc is a basis vector. Their actions onto 
the remaining basis vectors give zero. The corresponding matricially free annihilation 
operators are their adjoints pp^qiu). In some cases, it will be convenient to use the same 
notation even if hp^q{u) = 0, when we obtain trivial operators. 

Certain linear combinations of the matricially free creation and annihilation oper¬ 
ators are of special interest. We have stndied several matricial systems of operators 
constructed from the matricially free creation and annihilation operators or their sym¬ 
metrized counterparts when describing the asymptotic distributions of blocks or sym¬ 
metric blocks of Hermitian and non-Hermitian Gaussian random matrices [7,8]. We 
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define them below, mentioning also the corresponding ensembles of random blocks. We 
also give their realizations as operator-valued matrices, which identihes some of them 
with the generators of free group factors used by Voiculescu in [16]. 

(1) Non-trivial matricially free creation and annihilation operators can be realized 
as operator-valued matrices from the C^-algebra ^ (x) M„(C), where A is the 
C'*-algebra generated by the family {£(p,q,u) : 1 ^ p, g ^ r,u s Zl} of *-free 
creation operators living in the free Fock space. Namely, 


Pp,g(u) = £(p,q,u)0e(p,q) 

Pp,qiu) = ^{.PAaT ®e{q,p), 

where {e{p, g) : 1 < p, g < r} is the array of matrix units. If p is the vacuum 
state in this free Fock space and is the vector state on Afn(C) dehned by the 
basis vector of C”, then the state d'g can be identihed with p (x) 'ipq, as we 
showed in [11]. 

(2) Matricially free Gaussian operators are self-adjoint operators of the form 


^p,qiA^ Pp,qi^) 'P Pp,qiu) , 


and they are the canonical Gaussian operators in our framework. It turns out 
that they describe the limit distributions of unbalanced symmetric blocks of 
Hermitian Gaussian random matrices. 

(3) It is convenient to introduce the symmetrized creation operators as 

G = J + Pq,piA) if P < g 

[ Pq,q{u) if p = q 


and the symmetrized annihilation operators as their adjoints pp^q{u). 

(4) In order to describe the limit distributions of balanced symmetric blocks of 
Hermitian Gaussian random matrices, we need to use symmetrized Gaussian 
operators 

^p,qi'^) ~ Pp,qi'^) 'P Pp,qi'^) j 

where all operators are non-trivial. It is easy to see that if all creation and anni¬ 
hilation operators involved are non-trivial, then the above ones can be identihed 
with the operator-valued matrices of the form 


I dip, Q, u) ® e(p, q) + g{p, g, u)* (x) e(g, p) if p < q 
[ s(g,u) (x)e(g,g) if p = g 


where {g{p,q,u) : 1 ^ p ^ q ^ r,u s u} is a family of circular operators and 
{s{q, u) : 1 ^ q ^ r, u E Zl} is a family of semicircular operators. These matrices 
are generators of free group factors introduced by Voiculescu [16]. 

(5) In order to describe the limit distributions of the usual (non-symmetric) blocks 
of non-Hermitian Gaussian random matrices, we need to use matricial R-circular 
operators [11] 

Cp,q{'^) ~ Pp,qi'^ ) + Pq,p{u ) ) 

where the notation means that in order to construct one such operator labelled 
by u one needs two operators, which are labelled by u' and u", thus the index 
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set U has to be doubled. If the added operators are non-trivial, then we can 
identify the above ones with operator-valued matrices of the form 


Cp,qiu) = g{p,q,u)®e{p,q) 

for any 1 ^ p,q ^ r and any u e U, where {g{p,q,u) : 1 ^ p,q ^ r,u e Zl} 
is a family of circular operators (here, the covariances of pp,q{u) are symmetric 
and identical for u' and u"). Sums ({u) = YjpqCp,q{'^) are i?-cyclic matrices 
introduced in [13]. 

(6) In order to describe the limit joint distributions of the symmetric blocks of 
non-Hermitian Gaussian random matrices, we need to use matricial circular 
operators [8] which are natural symmetrizations of Cp^q{u), namely 

Vp^qi^') Pp,qi.^ ) Pp,qi.^ ) i 


where u e U and 1 < p < g < r. If the added operators are non-trivial, then 
the above ones can be identihed with the operator-valued matrices of the form 


I 9{p,q,u)®e{p,q) + g{q,p,u)®e{q,p) ii p < q 
\ s{q,u) ®e{q,q) li p = q 


where {g{p,q,u) : 1 < p, g < r, u e w} is a family of circular operators and 
{s{q, -u) : 1 < g < r, u e W} is a family of semicircular operators. 


3. Hermitian Gaussian Block Ensemble 

In free probability, we study the mixed moments of independent n x n random ma¬ 
trices {Y{u,n) : u E U} under the expectation of the normalized trace 

r(n) = E o Tr(n) 

where 

Tr(n)(yl) = -Tr(yl) 
n 

for any n x n matrix A as n oo. Let us recall the asymptotic freeness result of 
Voiculescu for independent Hermitian Gaussian random matrices [15]. 

Theorem 3.1. If we are given an ensemble of independent Hermitian n x n random 
matrices {Y{u,n) : u e Zl}, whose entries Yij{u,n) satisfy Yij{u,n) = Yj^i{u,n), are 
complex {0,1/n)-Gaussian if i A j and real {0,1/n)-Gaussian if i = j and the family 
{Yij{u,n) : 1 < f < j < n} is independent for any u, then 

lim Y{u, n) = uj{u), 

n—*’CC 

which should he understood as the convergence of mixed moments of matrices under 
T{n) to the mixed moments of the corresponding free standard semicircular operators 
u{u) under the vacuum state <h in the free Fock space. 

Remark 3.1. Let us make some remarks which will enable us to formulate our main 
results. 

(1) We stated the above theorem in a simplihed form, which will also be used when 
we state our results on the asymptotic distributions of random blocks. A more 
explicit formulation says that 

lim r(n)(E(ui,n). ..Y{um,n)) = . ..u{um)), 

n^co 



6 


R. LENCZEWSKI 


for any Ui,... ,Um ^ U, where {uj{u) : tt e w} is a standard free semicircular 
family, which means in particular that <h(a;(n)^) = 1 for any u. Alternatively, 
one could say that the family {Y{u, n) : u s U} is asymptotically free under 
r(n) and that the limit distribution of each Y{u, n) is the standard semicircular 
Wigner distribution W{0, 1) with density (27r)“^\/4 — on [—2,2], 

(2) In order to dehne a family of partial traces, take the decomposition of the set 
[n] := {1,2,..., n} into r disjoint intervals 

[n] := A^i u ... u Nr 

and denote by Uq the cardinality of Nq. Next, let 

I{n) = Di + ... + Dq 

be the corresponding decomposition of the n x n identity matrix /(n), that is 
{Dk)ij = 1 whenever i = j e Nk, with the remaining entries equal to zero, where 
1 < /c < r. Of course, the objects Nq,nq, Dq depend on n, but this is supressed 
in the notation. 

(3) By partial traces we then understand states of the form 

Tq{n) = EoTrg(n) 

where 

TTq{n){A) = -TliDqADq) 
riq 

and 1 < g < r. 

(4) By random blocks of a random matrix Y{u, n) we shall understand nxn matrices 
of the form 

Sp^q{u,n) = DpY{u,n)Dq, 

for any 1 < p, g < r and u e Zl, where r, n e N. In particular, if Y[u, n) is 
Hermitian, then Sq^p{u,n) = Sp^qiu^n)*. Clearly, we have the decomposition 

Y{u,n)= Yj 

l^p^q^r 

for any u, n. 

(5) We assume that all variables Yij{u, n) which belong to the same block Sp^q{u, n) 

have the same covariance E(hjj(M, n)) = Vp^q{u)/n whenever {i,j) e 

Np X Nq. We denote by V{u) = {vp^q{u)) the corresponding matrices of block 
covariances. Apart from the dimension matrix 

D = diag((ii, ...,dr), 

these matrices are additional parameters of the ensemble. In fact, the limit 
distributions will be expressed in terms of matrices B{u) = DV{u). 

(6) By random symmetric blocks of a random matrix Y{u, n) we shall understand 
nxn matrices of the form 

T - i + Sg^p{u,n) if p<q 

’ [ Sq^qlu,n) ifp = q 

for any 1 < p < g < r and u e Zl, where r, n e N. We have the decomposition 

Y(u,n)= Y Tp^qi'^,n) 

l^p^q^r 
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for any M,n. Clearly, if Y{u,n) is Hermitian, then Tp^g{u,n) = Tp^q{u,n)*. 

(7) When speaking of limit distributions of mixed moments of symmetric blocks, we 
will use a simplihed formulation, similar to that in the free case. For instance 

lim Tp^q{u,n) = ojp,q{u) 

in the Hermitian case will mean that 

rg(n)(Tpi,gi(Mi,n) .. .Tp^^g^{um,n)) 

= 'l/q(a;pj^qj (wi) . . ■ 0Jp^^q^{Um)) 

for any 1 < pi < gi < r,..., 1 < ^ 9 m ^ l, 1 < g < r and Mi, e Zl. 

The operators will always belong to one of the families dehned in Section 2. A 
similar formulation will be used for blocks Sp^g{u,n). 

We can describe limit joint distributions of blocks and symmetric blocks of Hermitian 
random matrices under rather general assumptions. It suffices to assume that the family 
{Y (m, n) : u e U} is asymptotically free and asymptotically free from {Di, ..., Dr} under 
r(n) and that their norms are uniformly boubded almost surely. This class includes 
unitarily invariant random matrices whose limit moments are compactly supported 
probability measures on the real line. This general version was proved in [8, Theorem 
6.1]. For the sake of simplicity, we restrict our attention here to the case of Hermitian 
Gaussian random matrices. 

Theorem 3.2. If we are given an ensemble {Y{u, n) : u s U} of independent Hermitian 
nxn random matrices whose entries Yij{u, n) satisfy Yij(u, n) = Yj^i{u, n), are complex 
{0,Vp^q{u)/n)-Gaussian if i ¥= j and {i,j) e Np x Ng, and real {0,Vg^g{u)/n)-Gaussian if 
i = j E Ng and the family {Yij{u, n) : 1 < i ^ j < n} is independent for any u, then 

lim Tp^g{u,n) = ojp^q{u) 

for any p ^ q and u e U, where convergence is in the sense of mixed moments under 
partial traces and the arrays {u)p^q{u)) are associated with the symmetric covariance 
matrices B{u) = DV{u), respectively. 

Sketch of the proof. First, let us assume that Vp^q{u) = 1 for any p, g, u. In that case, 
we can use asymptotic freeness of the family {Y{u, n) : u e Zl} and its asymptotic free¬ 
ness with respect to the family {Di,..., Dr} (these are deterministic diagonal matrices) 
to describe the limit joint distributions of the blocks Sp^q{u, n) under r(n) since 

Sp^q{u,n) = DpY{u,n)Dg 

for any p, g, u. Moreover, we know that the limit distribution of each Y[u, n) is the 
standard semicircle Wigner law IF(0,1) and a direct computation gives T{n){Dq) = 
Uq/n dq for any g. Now, it is not hard to show that the family of operators 

p,i 

where Up^q{u) = pp,q{u) + pp^giu)* for any u and pp^q{u) has covariance dp, is free with 
respect to T = '^gdg'^g. We abuse the notation a little since by u{u) we denoted a 
standard semicircular operator on the free Fock space, but this is justihed by the fact 
that each u{u) has the standard semicircle distribution under T. Therefore, we have in 
fact a realization of the standard free semicircular family {oj{u) : u e Zl} as operators 
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living in the matricially free Fock space A4. Now, it suffices to find an appropiate 
limit realization for (Di, ..., Df.}. We have shown in [8] that it is given by the family 
{Pi,..., Pr}, where 

Pq = l(^e{q,q), 

using the tensor product realization described in Section 2, since the family {u^u) : u e 
U] is free from {Pi,...,Pr} and d'(Pq) = dq. Therefore, any limit mixed *-moment 
of the random blocks Sp^q{u) under r(n) can be expressed as a mixed *-moment of 
the corresponding operators Ppu{u)Pg under T. Symbolically, since \min-^«)Y{u,n) = 
u){u), we have 

lim Sp^q{u,n) = Ppu{u)Pq = pp,q{u) + pq,p{uY := <,p,q{u) 

n^oo 

except that the moments of blocks are computed under r(n) and those of the operators 
Ppu{u)Pq are computed under T. In order to pass from r(n) to Tg{n), notice that 

Tl 

Tq{n){A) = — T{n){DqADq) 

Uq 

so the limit mixed *-moments of random blocks under Tq{n) are given by the mixed 
*-moments of the above operators under where 1 ^ q ^ r. All this holds provided 
dg A 0. The case when dg = 0 is slightly more complicated and is omitted here (we 
refer the reader to [8]). It easily follows that 

lim Tp^q{u,n) = ojp,q{u) 

for any 1 < p < g < r and u e U. It can also be seen that one can rescale blocks 
Sp,q{u,n), which means that one can rescale the block covariances and take arbitrary 
non-negative Vp^q{u) (except that we must have the symmetry Vp^q{u) = Vg^p(u) since 
Y(u,n) is Hermitian). This proves that in the case when the covariances are equal to 
Vp^q(u) within block Sp^g(u), the limit operator qp,g(u) gets rescaled by Vp^g(u) and thus 
the covariance of pp,g(u) becomes bp^g(u) = dpVp^g(u). This completes the proof. ■ 


Corollary 3.1. In particular, 

(1) if dp = 1 and dg = 0, then hm„^oo Tp^g(M, n) = Up^giu), 

(2) if dp = 0 and dg = 1, then lim^^oo Tp^g(M, n) = Uq^p{u), 

(3) if dp = 0 and dg = 0, then lim^^oo ^p,g('W, n) = 0. 

Proof. If dp = 1 and dg = 0, then Uq^p{u) = 0 since hq^p{u) = dqVq^p{u) = 0 and thus 

djp,q{u) reduces to Up^q{u). In turn, if Tp^q{u,n) is evanescent, then Up^q{u) = Uq^p{u) = 0 
and thus ujp^giu) = 0, which completes the proof. ■. 


Remark 3.2. It should be remarked that the matricially free Gaussian operators are 
not operatorial realizations of the usual (non-symmetric) blocks Sp^q{u) of Gaussian 
Hermitian matrices. In fact, it follows from the proof of Theorem 3.2 that 


lim Sp^q{u) 

n^oo 


g{p, q, u) (X) e(p, q) ii p ^ q 
s{q, u) (X) e{q, q) ii p = q 


for any 1 ^ p, g < r and u s U, where convergence is understood in the usual sense 
(mixed moments of blocks under partial traces Tg(n) converge to mixed moments of 
the corresponding operators under Tg). Here, {g{p,q,u) :l<p<g<r, wefijisa 
family of circular operators and {s{q,u) : 1 < g < r, m e W} is a family of semicircular 
operators and we assume that g{q,p,u) = g{p,q,u)* for p < q. 
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Example 3.1. Fix u e U (omitted in our notations) and p ¥= q. Let g{p,q) = ii + 
and thus g{q,p) = £* + £2 for some *-free creation operators £i,£ 2 - In view of Remark 
3.2, we get 

Ji^Tg(n)(S'g,p(n)S'p,g(n)S'q,p(n)S'p,g(n)) = ^ q{gq,pgp,qqq,pqp,q) 

= p{{£l + £2)i£i + £V){^i + ^2){£i + £ 2 )) 

= {dp + dpdq)Vp^q 

where p{£*i£i) = q{pl,qpp,q) = dpVp^q and ip{£l£2) = p{p*q,pPq,p) = dpVq^p and we use 
the symmetry Vp^q = Vq^p. 


4. Ginibre Block Ensemble 

Blocks of non-Hermitian Gaussian random matrices can be treated in a similar way. 
The associated ensmble of blocks cna be called the Ginibre Block Ensemble. The main 
difference is that the limit joint distributions of blocks are described by matricial R- 
circular operators Cp,q{u) and those of the symmetric blocks by matricial circular oper¬ 
ators Pp^q{u). 

Let us recall Voiculescu’s theorem on the asymptotic freeness of the ensemble of 
independent Gaussian random matrices with i.i.d. entries (Ginibre Ensemble) [15]. 

Theorem 4.1. If we are given an ensemble of independent n x n random matrices 
{Y{u,n) : u e U}, whose entries Yij{u,n) are complex {0,1/n)-Gaussian for any i,j 
and the family {Yij{u,n) : 1 ^ i, j ^ n} is independent, then 

lim Y{u, n) = g{u), 

n—*’O0 

where convergence is in the sense of mixed moments of matrices under T{n) to the mixed 
moments of the corresponding free standard circular operators g{u) under the vacuum 
state <h in the free Fock space. 

Let us now formulate an analogous theorem for blocks and symmetric blocks of 
Gaussian random matrices with i.b.i.d. entries. 

Theorem 4.2. If we are given an ensemble of independent n x n random matrices 
{Y{u,n) : u e u], whose entries Yij{u,n) are complex {0,Vp^q{u)/n)-Gaussian for any 
{i,j) E Np X Nq and the family {Yij{u,n) : 1 ^ i, j ^ n} is independent, then 

lim Sp^q{u,n) = Cp,q{u), 
n^oo 

where convergence is in the sense of mixed moments of blocks under partial traces Tq{n) 
to the mixed moments of the corresponding matricial systems of operators under the 
vacuum states Tg, respectively, in the matricially free Fock space. 

Sketch of the proof. The proof is similar to that of Theorem 3.1. In the case when all 
variables are i.i.d., we can use the asymptotic *-freeness of the family {Y{u, n) : u e Zl} 
under r(n) as n —>■ 00 as well as their asymptotic *-freeness from the family of diagonal 
matrices. One can check that the family {g{u) : u e U], where 

® 
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is *-free under 'I' as well as *-free with respect to {Pi,... ,Pr} under \h, where again 
Pq = 1 ® e{q, q). Moreover, each ri{u) has the standard circular distribution under 
'll = Xiq Since the asymptotic joint distribution of {Pi,..., Dr} under r(n) agrees 
with that of {Pi,..., P^} under \h, we must have 

lim Sp^q{u,n) = lim DpY{u,n)Dq = Ppr]{u)Pq = Cp,q{u) 

n—*’O0 n—>(X) 

where convergence is understood as described in Remark 3.1, which completes the proof. 


Corollary 4.1. Under the assumptions of Theorem 4-2, it holds that 

lim Tp^q{u, n) = r]p,q{u) 
n^oo 


for any p, q, u 


Proof. This is an easy consequence of Theorem 4.2. 


Example 4.1. Computations of limit mixed (*-) moments of blocks reduce to the com¬ 
putation of moments involving matricial P-circular systems of operators. For instance, 

lim Tq{n){Sp^q{nySq^p{n)*Sq^p{n)Sp^q{n)) = ^qiCp,qClpCq,pCp,q) 


dpdqVp^qVq^p 


where p ¥= q, since 


Cp,q = (^1+ '^2)®e(p,g), Cg,p= {^3 + £l)®e{q,p) 
where {£i, ^ 2 , ^ 3 i ^ 4 } is a *-free system of free creation operators with covariances ^^(^i^i) = 
dpVp^q and = dqVq^p. In fact, it is well known that any pair of free circular oper¬ 
ators can be written in the form £1 -I- ^3 + respectively. 


5. Products of independent Gaussian random matrices 


The first concrete application of our method concerns products of independent rect¬ 
angular Gaussian random matrices [10]. For any given p e N and any n e M, consider 
the product of independent rectangnlar Gaussian random matrices 

P(n) =Xi(n)X 2 (n)...Xp(n), 

where n e N and all entries of these matrices are assnmed to be independent ( 0 , 1 /n)- 
Gaussian variables. If Xj{n) is an Uj-i x Uj matrix for any 1 < j < p, we assnme 
that 

lim — = dj 

for any j e {0,... ,p} (it is convenient to start with uq rather than with ni). Let ro(n) 
be the trace over the set of hrst no basis vectors composed with classical expectation. 

Theorem 5.1. Under the above assumptions, it holds that 

j^ro(n) (^(5(n)5*(n))^j = Pk{do,di, ... ,dp), 


Pki^do^ di^ • • •; dpj ^ j 

jo+...+jp=pk+l 


1 fk 
k\jo 


Jp 




where 
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for any natural k and jo, ji, ■ ■ ■, jp- These polynomials are called multivariate Fuss- 
Narayana polynomials. 

Sketch of the proof. The proof is based on embedding the matrices Xi{n),..., Xp{n) 
in symmetric blocks Ti 2 ( 11 ),..., Tp p+i(n), respectively, of a large square Gaussian ran¬ 
dom matrix Y{n) of dimension NxN, where N = ni + .. .+np+i, with (0, l/n)-Gaussian 
entries for any n (we use only one matrix and thus we omit u in our notations). Gom- 
puting the moments of B{n)B*{n) under ro(n) becomes now the normalized partial 
trace over the first subset of basis vectors. We then use Gorollary 4.1 to realize the 
limit moments in terms of 771 ^ 2 , • • • ,'hp,p+i and their adjoints. These limit moments can 
be computed explicitly, which was done in [10], which completes the proof. ■ 


Let us make some additional remarks on the above result: 


(1) The special case of p = 1 corresponds to Wishart matrices. If we set do = I and 
di = t, we obtain 


Pk{l,t) 





the moments of the Marchenko-Pastur distribution with shape parameter f > 0, 
namely 

'J {x — a){h — x) 

Qt = maxjl - f, Ojdo H- - \\a,b\[x)dx 

ZTTX 

where a = (1 — and 6 = (1 -I- 

(2) The moments of the Marchenko-Pastur distributions are known to have the form 
of Narayana polynomials 



for any k s N. It is easy to show that Pk{l,t) = Nk(t), but our formula is more 
suitable for multivariate generalizations. 

(3) If do = 1, the multivariate Fuss-Narayana polynomials become the moments of 
the p-fold free multiplicative convolution of the Marchenko-Pastur distributions 
with arbitrary shape parameters. Namely, 


Pfc(l,ti, ...,tp) = rukigti 0 ptj 0 • • • 0 Qt^) 

where mk{fi) stands for the A:th moment of p and 0 stands for the free multi¬ 
plicative convolution. 

(4) A less general class of polynomials was studied in the context of free Bessel laws 

TTpp = 0 Qt, 

where p e N, dehned by Banica, Belinschi, Gapitaine and Gollins [2]. The 
moments of free Bessel laws are polynomials in one variable t, 
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called Fuss-Narayana polynomials. Our polynomials are natural multivariate 
generalizations of these polynomials. Clearly, Qk{t) = Pk{l,... where 1 

appears p times. 


6. New random matrix models 


The results on the asymptotics of random blocks can be applied to the construction 
of some new random matrix models, as we showed in [9]. In this paper, we choose 
to present a simple version of such a model for monotone independence [12]. In a 
similar way, one can construct random matrix models for boolean independence an 
s-free independence (freeness with subordination). In fact, more general versions, not 
restricted to Gaussian random matrices have been constructed in [8]. 

Monotone independence will appear in the study of the asymptotic joint distributions 
of two independent Hermitian Gaussian random matrices of the same block form. Our 
assumptions are the following: 

(Al) We have a family of independent Gaussian random matrices 

, N f A{u,n) B{u,n) \ 

^ ("■")=(cKn) D{u,n) ) 


where u e Zl, and blocks {A{u,n)) are evanescent, symmetric blocks built from 
{B{u,n)) and (C'(u,n)) are unbalanced, and blocks {D{u,n)) are balanced, 
(A2) the matrices are Hermitian, thus the off-diagonal blocks are Hermitian conjugate 
and the diagonal blocks are Hermitian 

(A3) the complex Gaussian variables hjj(u,n) have zero mean and have identical 
covariances within blocks, namely 

E{Yij{u,n)Yij{u,n)) = 

whenever the pair {i,j) belongs to the block indexed by (p, g), 

(A4) the decomposition of the identity matrix corresponding to the block decompo¬ 
sition is given by 

I (A) = Di + D2 

for any n e N. 


Theorem 6.1. Under assumptions (Al)-(A4), if Zl = {1,2} and Vp^q{u) = 1 for any 
p, q, u, the pair {B{1, n)-\-C{l, n), Y (2, n)} is asymptotieally monotone independent with 
respeet to the partial trace Ti{n). 


Proof. By Theorem 3.2, the proof reduces to showing that the pair 

{ 1 ^ 2 , 1 ( 1 ), <^ 2 , 1 ( 2 ) + 012 , 2 ( 2 )} 

is monotone independent with respect to since, by assumption, the asymptotic 
dimensions are di = 0 and ^2 = 1, which means that the remaining operators can be 
neglected. Denote a = 022 , 1 ( 1 ) and b = 022 , 1 ( 2 ) -1- 022 , 2 ( 2 ). We need to show that 

i{wiaibia2W2) = '^i{bi)'^i{wiaia2W2) 

for any 01,02 e C[o, li], bi e C[6, 12 ], where li = l 2 ,i and I 2 = 1 and wi,W 2 are 
arbitrary elements of C(a, 6, li, I 2 ). It suffices to consider the action of a and b onto 
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their invariant subspace in A4 of the form 


M' = cffi © (j-(2) (g) n(i)) © (j-(2) © n(2)) 


where J^(2) = ^^(© 62 , 2 ( 2 )) with the vacuum vector hi and ?{(u) = © 62 , 1 ( 11 ) for m e {1,2}, 
where we identify © 62 , 1 ( 11 ) with 62 , 1 ( 11 ). Now, the range of any polynomial in a is 
©ffi © 'H(l) since 


J hli if k is even 
I 62 , 1 ( 1 ) if k is odd 


and lihli = hli, 1 i 62 ,i( 1 ) = 62 , 1 ( 1 ). Therefore, it suffices to compute the action of any 
polynomial in b onto hli and 62 , 1 ( 1 ). Now, the action of powers of b onto hli and onto 
62 , 1 ( 1 ) is the same as the action of the free Gaussian operator onto the vacuum vector 
in the free Fock space. Namely, we have 

b^^Qi = Ckkli mod {M'eCQi) 


and 

62 " 62 ,i( 1 ) = ©^ 62 , 1 ( 1 ) mod (M'©(©Gi©?^(l)) 
where Ck is the kth Catalan number and 


= 62^-^e2,i(l) = 0 mod (M' © (©Gi © 7^(1)) 

for any k s N. Thus Ti(6^^) = Ck and, moreover, since M’ © (©Gi © 'H(l)) c Kera, 
the required condition for monotone independence holds if bi is a positive power of b. 
It is easy to see that it also holds if bi = I 2 , which completes the proof. ■ 


Interestingly enough, unbalanced symmetric blocks were also used in the construction 
of a simple random matrix model for free Meixner laws [1,4]. These are probability 
measures on the real line associated with the sequences of Jacobi parameters of the 
form 

(Q!i,Q! 2 ,a 2 ,.. •) and (/5i,/ 32 , ^ 2 ,...) 

and thus we can asy that they are associated with quadruples (ai, a 2 ,/3i,/ 32 ). Let us 
formulate the theorem in the most interesting case when /3i and (32 are positive. 

Theorem 6.2. Under assumptions (A1)-(A4), let (3i{u) = U 2 ,i(m) > 0 and (32{u) = 
V 2,2 > 0, for any u e Zl. Then 

(1) the asymptotic distributions of the matrices 

M{u, n) = Y{u, n) + ai{u)Di{n) + a 2 {u)D 2 {n) 

under the partial trace Ti{n) are the free Meixner distributions associated with 
the parameters (q:i(m), 0 : 2 ( 11 ), /5i(m), / 92 (ii)); respectively, 

(2) the family {M{u, n) : u e U} is asymptotically conditionally free with respect to 
the pair of partial traces (ri(n), T 2 (n)) as n ^ co. 

Sketch of the proof. It follows from Theorem 3.2 that 

lim Y{u, n) = 002 ,i{u) + uJ 2 , 2 {u) 

where the variances of U 2 ,i{u) and U 2 , 2 {u) are 62 . 1 ( 11 ) = 1 ^ 2112 , 1 ( 11 ) = /3i{u) and 62 , 2 ( 11 ) = 
d 2 V 2 , 2 {u) = /32 {u) since di = 0 and d 2 = 1- The fact that di = 0 is the reason why 
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and 011^2(w) become trivial operators and that is why they do not show np in 
the limit realization. Therefore, 

lim M{u,n) = 7(11) := uj2,i{u) + a;2,2(n) + ai{u)Pi + a2{u)P2, 

n-i-co 

where Pi and P2 are as in the proof of Theorem 3 . 2 . This implies, in particnlar, that the 
asymptotic distribntion of M{u,n) nnder Ti(n) agrees with that of 7(11) nnder Ti. The 
proof that the moments of 7(11) are the moments of the free Meixner law associated with 
the parameters (0:1(11), 02(n),/ 5 i('u),/32('u)) is pnrely combinatorial and can be fonnd in 
[ 9 ]. We also refer the reader to [ 9 ] for the proof of asymptotic conditional independence 
of the family of snch matrices. ■. 


References 

[ 1 ] M. Anshelevich, Free martingale polynomials, J. Fund. Anal. 201(2003), 228-261. 

[2] T. Banica, S.T. Belinschi, M. Capitaine, B. Collins, Free Bessel laws, Canad. J. Math. 63 (2011), 
3-37. 

[3] F. Benaych-Georges, Rectangular random matrices, related convolution, Probab. Theory Relat. 
Fields 144 (2009), 471-515. 

[4] M. Bozejko, W. Bryc, On a class of free Levy laws related to a regression problem, J. Fund. Anal. 
236 (2006), 59-77. 

[5] J. Cuntz, W. Krieger, A class of C*-algebras and topological Markov chains. Invent. Math. 56 
(1980), 251-268. 

[6] R. Lenczewski, Matricially free random variables, J. Funct. Anal. 258 (2010), 4075-4121. 

[71 R. Lenczewski, Asymptotic properties of random matrices and pseudomatrices, Adv. Math. 228 
(2011), 2403-2440. 

[8] R. Lenczewski, Limit distributions of random matrices, Adv. Math. 263 (2014), 253-320. 

[9] R. Lenczewski, Random matrix model for free Meixner laws, Int. Math. Res. Notices (2014), 
rnu041, 26 pages. 

[10] R. Lenczewski, R. Salapata, Multivariate Fuss-Narayana polynomials with appplication to random 
matrices, Electron. J. Combin. 20 (2013), Issue 2 (electronic). 

[11] R. Lenczewski, Matricial i?-circular systems and random matrices, arXiv:1311.6420 (2014). 

[12] N. Muraki, Monotonic independence, monotonic central limit theorem and monotonic law of small 
numbers, Infin. Dimens. Anal. Quantum Probab. Relat. Top. 4 (2001), 39-58. 

[13] A. Nica, D. Shlyakhtenko, R. Speicher, i?-cyclic families of matrices in free probability, J. Funct. 
Anal. 188 (2002), 227-271. 

[14] D. Shlyakhtenko, Random Gaussian band matrices and freeness with amalgamation, Int. Math. 
Res. Notices 20 (1996), 1013-1025. 

[15] D. Voiculescu, Limit laws for random matrices and free products. Invent. Math. 104 (1991), 
201 - 220 . 

[16] D. Voiculescu, Circular and semicircular systems and free product factors. Progress in Math. 92, 
Birkhauser, 1990. 

Romuald Lenczewski, 

Instytut Matematyki I Informatyki, Politechnika Wroclawska, 

Wybrzeze Wyspianskiego 27, 50-370 Wroclaw, Poland 

E-mail address: Romuald.Lenczewski@pwr.wroc.pl 


